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Constraints on the decay of dark matter to dark energy from weak 
lensing bispectrum tomography 
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ABSTRACT 

We consider a phenomenological model for a coupling between the dark matter and dark en- 
ergy fluids and investigate the sensitivity of a weak lensing measurement for constraining the 
size of this coupling term. Physically, the functional form of the coupling term in our model 
describes the decay of dark matter into dark energy. We present forecasts for tomographic 
measurements of the weak shear bispectrum for the DUNE experiment in a Fisher-matrix 
formalism, where we describe the nonlinearities in structure formation by hyper-extended 
perturbation theory. Physically, CDM decay tends to increase the growth rate of density per- 
turbations due to higher values for the CDM density at early times, and amplifies the lensing 
signal because of stronger fluctuations in the gravitational potential. We focus on degenera- 
cies between the dark energy equation of state properties and the CDM decay constant rele- 
vant for structure formation and weak lensing. A typical lower bound on the CDM decay time 
^ 1.1 /Ho = 75.3 Gyr/Zi which could be provided by DUNE would imply that it would not 
possible to produce the dark energy content of the universe by CDM decay within the age of 
the Universe for a constant equation of state parameter of w close to - 1 . 
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1 INTRODUCTION 

Dark energy as a cosmological fluid is evoked for explain- 
ing the late-time cosmic acceleration, which has been observed 
in a number of channels, e.g. the cosmic microwave back- 



ground (CMB) anis otropies jSpergel et al. 2003 ) and the integrated 
Sachs-Wolfe-ef fect (Boughn & Crittende r]|200a:lNolta et aLlllOOSl : 
iGiannantonio e t al. 2006, 2008; Rassat et al 2007h . The numerical 
similarity of the values of the dark energy density Q.^ ^ 0.75 and 
the matter density Q.„, ^ 0.25 constitues the coincidence problem: 
Why is the ratio Q.^/Q.,„ = 0(1) today, or, why is the Hubble ex- 
pansion dominated by dark energy after the formation of galaxies 
and the large-scale structure? 

We focus on coupled models of dark energy proposed by (see 
IValiviit j[2003 : iBoehmer et al.ll2008l) . where the dark energy (DE) 
is generated by decay of dark matter (CDM). The CDM decay rate 
r is a constant free parameter, in addition to the matter density n„ 
and the equation of state w of the dark energy. The time evolution 
of the matter density p„, and the d ark energy density is described 
by the energy balance equations teoehmer et al.ll2008h . which can 
be written in matrix form as 



3H + r 
-r 





3H(l +w) 



Pm 

P<l> 



0, 



(1) 
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with the Friedmann constraint = | (p,„ +p^)- Note that the 
baryons are not coupled to dark energy, since this would be sub- 
ject to sever constraints from 'fifth-force' experiments. We will ne- 
glect the baryons (and the radiation) in our analysis, which should 
not be a serious limitation, however, due to the fact that the baryon 
fraction /j = Oi/Q,,, ==0.16 has a small numerical value. 

We also assume a homogeneous dark energy fluid, even 
though the decay of CDM in structures will naturally lead to local 
overdensities of dark energy. If the dark energy sound speed has a 
value close to 1 , the dark energy will diffuse rapidly away from the 
dark matter structures and constitute a homogeneous fluid. Further- 
more, the density fluctuations measured by weak lensing are only 
mildly overdense such that this assumption may be justified. Fur- 
thermore, we will work with the Newtonian Poisson equation an d 
not the full general relativistic expression ( l01ivares"et ai] |2006h . 
which is justified as the dark energy is assumed to be homogeneous. 

iBoehmer et akl ( |2008|) provide a thorough discussion of the 
background dynamics of such a cosmology in the case where the 
dark energy is an exponential quintessence field. In order to anal- 
yse the growth of structure and weak lensing, we will use a sim- 
plified, phenomenological model of dark energy by introducing a 
parametrised time-variable equation of state. We aim to model a to- 
mographic weak lensing measurement, and thus provide forecasts 
for parameter constraints and to quantify parameter degeneracies 
as one can expect from a deep weak lensing survey. 

Gravitational lensing, in particular, has been shown to be 
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a very powerful observational probe for investigating the in- 
fluence of dark energy on structure formation and the ge- 
ometry of the univ erse ( Schneider et al.' 'l992; Meliie] 1 19991 : 
iBartelmann & Schneid er 2001; Refregier 2003). even in the 
nonlinear regime of s t ructure formation jjain & Seliad 1 19971 : 
ISernardeau et alj 1 19971 : iBenabed & Bemardeaij l200lh . Lensing 
data is best used in tomographic measureme nts for constrain- 
ing dark energy equation of state properties JHuI 1 19991 |2002| : 
iHeavenj l2003l : Uain & Tavloil [2003) . where one either measures 
the power spectrum or the bispectrum of a weak lensing quan- 
tity j Kilbinger & Schneidei"2005'; 'Schneider & Bartelmamj|l997l : 
iBern stein & Jain 2004; Dodelson & Zhana 20o3)- Supplementing 
the recent paper bv lLa Vac ca & Colombo ( 2008), who derived lens- 
ing bounds on interacting models from weak lensing power spectra, 
we focus on bispectrum tomography, and we use a more general, 
albeit phenomenological cosmological model. Bispectra have the 
advantage that the perturbative treatment is easier to carry out and 
that they are sensitive on the transition from linear to nonlinear dy- 
namics in structure formation. 

After introducing the cosmological model and the peculiarities 
of gravitational lensing in the decaying CDM models in Sect. [J] 
we compute the weak lensing bispectrum and tomographic mea- 
surements in Sect. [3] Fisher-constraints on cosmological parame- 
ters are derived in Sect. |4] and the main results are summarised in 
Sect. |5] The parameter accuracies are forecast for the weak lens- 
ing survey proposed for the Dark UNiverse Explore^ (DUNE). For 
the fiducial model, we take a spatially flat ACDM cosmology with 
w = -1, adiabatic initial conditions and stable CDM (F = 0). Spe- 
cific parameter choices are Hq = lOO/i km/i/Mpc with h = 0.72, 
n„ = 0.25, CTg = 0.8 and «j = 1. 



2 GROWTH FUNCTION AND WEAK LENSING 
2.1 Decaying dark matter 

The decay constant F will be expressed in units of the Hub- 
ble constant Hq. We rewrite eqn. ([TJ with the scale factor 
a as indep endent variable, and in t roduce the common param- 
eterisation ^Turner & Whit3 1 1991 IChevallier & Polarskil 1200 ll : 
iLinder & Jenkinsll2003h . 



w(a) = Wo -I- (1 - a)Wa. (2) 

The Hubble function H(a) is given by solving the differential equa- 
tion with solutions for p„,(a) and p^(a), 

= - {pm(a) + p^(a) [1 + w{a)]) . (3) 

The initial conditions are set at the present epoch ao = 1, and we 
integrate backwards to a = 10"-. With the solution of H{a), the 
comoving distance is 



1 a^H(ay 



(4) 



and the density parameters can be defined by normalising 0„(a) = 
Pm(a)/Pcrii(a) and n^(a) = p^(a)/pc^i,(a) with the critical density 
Peril = 3H^{a)/(8nG). Table[T]gives an overview of the dark energy 
models considered in this paper. 

The Hubble function and its first derivative are depicted in 
Fig.[T] where the SCDM-scaling is divided out. 
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Table 1. Summai'y of the four primai'y dai'k energy models considered in 
this paper, where F is the CDM decay constant. 
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Figure 1. The scaled Hubble function Ha(a) = a ' H{a)/HQ (thick lines) 
and the derivative H\(a) = a^^^dH(a)/da/Ho (thin lines) in models with 
decaying CDM, in comparison to models with stable dark matter: ACDM 
(solid line), ApCDM with F = | (dashed line), i^CDM (dash-dotted line), 
and (^rCDM with V = \ (dotted Une). 



d" H(a) 



da" Ho " 

An interesting feature, which will impact on the solution of the 
growth equation and on the Poisson equation, is the faster scaling 
of the Hubble function with the scale factor a in models with de- 
cay. This arises because the CDM density decreases not only by 
redshifting, which would be oc a^^, but also by decay. In addition, 
CDM decay causes the Hubble function to vary more gradually, 
which similarly can be achieved by introducing a variable equa- 
tion of state of the dark energy fluid. Looking at the scaled first 
derivative of the Hubble function Hi{a) illustrates that the inter- 
play between CDM decay and a variable dark energy equation of 
state gives rise to qualitatively new features in the Hubble function 
at scale factors 0.2 < a < 0.5 , compared to dark energy models 
with stable CDM. 



2.2 Structure formation 

The linear growth function D+(a), describing the homogeneous 
growth of structure, 6(x,a) = D+{a)S(x,l) under Newto- 
nian gravity, is obtained by solving the diffe rential equation 
jWang & Steinhardlll99i : lLinder & Jenkinsll2003h , 



d^ 1 / din// 

—£»+(«) + - 3 + -— 
aa-^ a\ a ma 



^DAa) = -^n,„(a)/)+(a), 
da 2a^ 



(6) 



which is still valid in decaying CDM models, as the decay does not 
affect the overdensity field 6(x) = (p(x) - (p))l(p), (p) = Cl„,Pciif 
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Figure 2. Growth function D+{a) for four dark energy models: ACDM 
(solid Une), ArCDM with r = 5 (dashed hne), i/iCDM (dash-dotted line), 
and (ArCDM with T = ^ (dotted line). Additionally, the growth function 
D+(a) = a for SCDM is plotted (solid straight line). 
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Figure 3. Lensing efficiency functions W{x)lx for the four dark energy 



models: ACDM (solid line), ApCDM with T 



(dashed line), 0CDM 



(dash-dotted line), and ^pCDM with T ■ 
sion into tomography bins. 



(dotted line), without subdivi- 



The initial conditions for decaying CDM models are different com- 
pared to those of standard dark energy models, as the scaling of the 
Hubble function as well as the time-evolution of Q.,„ are changed 
by the decay. The asymptotic behaviour of the growth equation can 
be separated out by assuming oc a" with a positive constant 

a at early times. Substitution into eqn. {SJ yields a quadratic equa- 
tion for a, which can be solved for an expression of a depend- 
ing on 3 + dln/f/dlna and n„(a), evaluated at the initial time 
a\. From that, one obtains the initial conditions £)+(«]) = a" and 
dD+/da = orap'. 

Fig.|2]shows the growth functions D+(a) (normalised to unity 
today) in the four cosmologies considered. Evolving dark energy 
has the property of suppressing structure formation at an earlier 
time, which can be partially compensated by CDM decay (be- 
cause the gravitational fields generated by the overdensity 6{x) are 
stronger if Q,„(a) has a higher value), hinting at degeneracies be- 
tween the equation of state parameters and the CDM decay rate. 
The ACDM and 0rCDM-models, for example, are almost indis- 
tinguishable. The influence of the two terms 3 d In HjA In a and 
n„,(a) on the evolution of the growth factor D+{d) is discussed in 
Appendix lAl 



Inclusion of the redshift distribution of background galaxies and 
reformulating the integration yields: 



where we abbreviated 
G(x)= r dx' P(z)-P;^'~^ 



(9) 



(10) 



d^' X' 

with the redshift distribution p{z)Az of the lensing galaxies. From 
that o ne finds an an alogous expression for the convergence spec- 
trum ( lLimbej|l954t) . 

9 rxH 

CM = ^ ixG(xf[aH(a)fnl(j^)Dl(a)P(k = f/x), (U) 

which results from projection of the CDM spectrum P{k). The stan- 
dard results are recovered by setting 



(12) 



which is valid only for stable CDM-models. The lensing efficiency 
function W(x) can be isolated from the above expresssions, yield- 



2.3 Weak lensing 

The non-standard scaling of the CDM density with scale factor 
makes it impossible to carry out a number of simplifications when 
deriving the weak lensing convergence. Substituting the critical 
density Pc.,it(a) = 3H^(a)/(SnG) with (p){a) = Q.,„(a)pcnt(a) into the 
expression AO = 4nGa^ {p){a) S for the comoving Poisson equation 
yields 



AO : 



-a^HHa)il,„{a)6, 



with explicit 
pression one 



functions H{a) and Q.„ 
obtains for the weak 



( iBartelmann & Schn eider 2001) 



Jo 



dx' a^H^(a)Q.„(a) 



XX -X' 



X 



(7) 

(a). Using this ex- 
lensing convergence 



(8) 



mx) 



2c2 



a^H\a) a„,{a) G{x)x. 



(13) 



Lensing efficiency functions Wix)lx for the four exemplary cos- 
mologies are given in Fig. [3] At low redshifts, models with evolv- 
ing dark energy attain higher values for Wix) compared to models 
with a cosmological constant, and at high redshifts, models with 
decaying CDM have higher values for Wix) compared to mod- 
els with stable CDM, with an interesting crossing at a distance of 
^ 100 Mpc/A. This behaviour is caused by the evolution of H{a) 
as discussed in the previous section and indicates that it is possible 
to increase the redshift at which the lensing signal originates by in- 
creasing r and by choosing an equation of state model with mean 
w close to -1. 

Fig.|4]illustrates the impact of CDM decay on linear conver- 
gence power spectra: Models with decaying CDM exhibit larger 
values for the power spectrum, which is due to the fact that in these 
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Figure 4. Linear convergence power spectra Cuit) for the four specified 
dark energy models: ACDM (solid line), ApCDM with F = i (dashed Hne), 
(^CDM (dash-dotted line), and ^pCDM with T = \ (dotted line), without 
subdivision into tomography bins. 



models the matter density in the past was higher compared to mod- 
els with stable CDM, which leads to stronger gravitational poten- 
tials and hence a stronger light deflection, which is described by the 
higher values of the lensing efficiency function W(x), with a minor 
effect from the growth function £)+(«). This feature might in fact 
be able to reconcile the high value for erg required by weak lensing 
surveys with the erg value following from CMB data, by generating 
high values for the lensing signal with a comparatively small value 
for CTg, if one allows for CDM-decay. 

In summary, the lensing signal measures a combination of 
structure growth D+(a), coupling strength n^(a) and geometrical 
factor j^'Cfl), which all are influenced by the interaction term in the 
evolution equations for the dark matter and dark energy density 
and make up the form of the lensing efficiency W(x) and the finally 
the expression for the spectra C^(t) and 6,(f|, ^s)- The epoch- 
varying coupling of the light to the gravitational potential, due to 
the non-standard Poisson equation e qn. (|7) is an important mecha- 
nism, which is not clearly treated in lLa Vacca & Colombd ( 1200 Sl) . 
It is worth noting that models with interacting dark fluids constitute 
a new class of models in this respect, in comparison to models with 
non-interacting evolving dark energy. 

Concerning the numerical evaluation, we generate a lookup ta- 
ble with 200 values of the relation between scale factor and comov- 
ing distance a(x), the Hubble function H(a), the growth function 
£>+(«), the matter density parameter Q.^i'^), and finally the nonlin- 
ear wave vector k„\i„ for a given cosmological model and use cubic 
splines for interpolating between these values. Additionally, values 
for the tomography lensing efficiency functions G,Ct) for each red- 
shift bin i are cached, as a function of comoving distance. These 
measures provide a significant speed-up in all computations, be- 
cause we do not have to reevaluate the densities p,„ and for each 
time step, to solve for the Hubble function H(a) and the density 
parameter n,„(a), and use these results for computing the growth 
function D+(a), the comoving distance and the lensing efficiency 
function Gi(x)- 



3 WEAK LENSING BISPECTRUM TOMOGRAPHY 

In this section we compile the necessary formulae for a tomo- 
graphic measurement of the bispectrum and for a Fisher-matrix 
analysis in order to derive constraints on cosmological parameters 
including the CDM-decay rate T and to quantify parameter degen- 
eracies, especially between T and the dark energy properties wq and 
Wa, from a weak lensing observation. 

We choose to consider a measurement of the bispectrum 
fi)t(^i, ^2, ^s) rather than a measurement of the spectrum C;t(^) be- 
cause both yield very similar constrain t s on d ark energy parame- 
ters, as shown bv lTakada & JaiiJ l l2003bl , l2004) . have a high covari- 
ance (which is difficult to quantify as it involves the computation 
of a 5-point correlation function) and are by no means independent 
measurements, and because of the fact that perturbation theory for 
describing the nonlinear evolution of the cosmic density field is in 
fact simpler for the 3-point functions compared to the 2-point func- 
tions: A consistent description of the nonlinear corrections to the 2- 
point function would involve contributions from both convergence 
fields perturbed to first order, and from one field perturbed to sec- 
ond order while the other field remains unperturbed. Contrarily, the 
leading contribution for the bispectrum involves the perturbation 
of a single field to first order, and decomposition of the resulting 
4-point function with the Wick-theorem into a product of 2-point 
functions. 

The most important difference between the spectra and bis- 
pectra concerning the line-of-sight integrations, however, is the fact 
that the spectrum is proportional to W^(x)D-l(x), whereas the bis- 
pectrum measures different powers of the lensing efficiency func- 
tion and the growth rate, namely W^(x)D'l(x)- 



3.1 Perturbation theory 

For the linear power spectrum P(k), which describes the fluctua- 
tion amplitude of the density field 6(k), {6(k)6(k')) = (Irrfdoik + 
k')P(k), we make the ansatz 

P(k)cc k"'T~(k). (14) 

The shape of the transfer fu nction T(k) is well app roximated by the 
fitting formula suggested bv lBardeen et alj|l986l) . 

ln(l + 2.34g) 



23Aq 



■ (l + 3.89^ + (16.1g)- + {5A(>qf + {6.1\qf) 



wh ere the wave vect or k = qF, is rescaled with the shape parameter 
T, dSugivamall 19951) , 



: Q.„h exp 



-nJ 1 + 



(15) 



The fluctuation amplitude is normalised to the value erg on the scale 
R = 8 Mpc/h, 



4 = ^ J dkk'W\kR)P(k), 



(16) 



with a Fourier-transformed spherical top-hat W(x) = 3j](x)/x as 
the filter function. jr(x) denotes the spherical Bessel function of 
the first kind of order £ ( Abramowitz & Stegun 1972). 

The first order contribution to the bispectrum Bs(ki,k2, k ^) of 
the den si ty field from nonlinear structure formation is given by l lFrvl 
ll984allbllT^ada & Jainll2003ah : 



Bs{kt,k2,k3)= M(k„kj)PML(kdPNL(kjh 



(17) 



(>J)6|1.2.3| 
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with the classical mode coupling functions, 



10 Iki kj\ 4 



4 , 



(18) 



which is replaced by the mode coupling function in hyper-extended 
perturbation theory, 



10 



k, , kj 



M{k„kj) = —a{k,)a{k,)+b{k,)b{k,)\j + -ijx+-c(^,)c(^X.(19) 

where x = kikjl [kik^ denotes the cosine of the angle between 
kj and kj. The coefficients a(k), b(k) and c(k) are given b y 
dScoccimarro & Frieman|[T99^: ls coccimarro & Couchmar]l200lh : 

1 + (T^°Hz) V0.72(m) (^/4)"+3-5 



a(k) 



b(k) 



c(k) 



1 + (g/4)'"+3-5 

1 + 0.4(n + 3)<7"+^ 



1 + 



1 + 



4.5 



(2q)" 



(20) 



(21) 



(22) 



1 + (2g)''+3-5 

where the time evolution of the fluctuation amplitude is given by 
the linear growth function, crgCz) = D+(z)o-s. The wave vectors are 
rescaled with the nonlinear wave number, q = kikm,. The logarith- 
mic slope of the linear power spectrum. 



n{k) 



A\T\P(k) 



(23) 



dln/t 

can be straightforwardly determined with the polynomial fit to the 
transfer function T{k). The nonlinear wave number at scale factor 
a is given by the scale at which the variance cr of the density fluc- 
tuations becomes unity, 

a- = I d?kDl{z)P{k) = 1 4nki, Dl(z)P(kt,L) = 1. (24) 



d'kDl(z)P(k)=l-^4nkli^Dl( 





Finally, the saturation parameter Q(n) can be computed from the 
logarithmic slope n of the linear CDM spectrum. 



Q(n) 



4-2" 



(25) 



1 + 2'"+' ■ 

Furthe rmore, we use the parameterisation proposed bv lSmith et alj 
( I2OO3 ) for the nonlinear CDM spectrum -Pnl(^) and its slow time 
evolution, which is particularly useful as the time evolution is pa- 
rameterised with Q.,„(a). 



3.2 Angular bispectra 

The spherical bispectrum Bi^((i, ti, is related to the flat-sky bis- 
pectrum B,{ti J2,ti)^y CMiralda-Escude, 1991 : .Kaiser. 1992) 











WU^'it,, + 1) 

4;r 



S,(f,,f2,^3),(26) 



where 

/A 








AxPcAx)Phi^)Phi^h 



(27) 



X = cos Q, denotes the Wigner-3 j symbol, which results from in- 
tegrating over three Legendre polynomials Pt(x). The Wigner-3_/ 
symbol cancels non-admissible configuratio ns which would violate 
the tr iangle inequality - t}^ < 4 < ^/-H^,- 1 Abramowitz & Stegurj 
Il972h . We use the Limber-equation (lLimbeilll954l) in the flat-sky 
approximation. 
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Figure 5. Equilateral convergence bispectra Bg(f): ACDM (solid line), 
ArCDM with T = ^ (dashed line), <pCDM (dash-dotted line), and 0rCDM 
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W\x)D\{x)BAkuk2,k,), 



(28) 



with kp = tplx, p = 1,2,3, for projection of the angular con- 
vergence bispectrum B^. The factorials in the Wigner-3 j symbol 
are evaluated using the Stirling-approximation for the F-function, 
F(n + 1) = n! with 



F(x) ^ V27r exp(-x) r'"i 



(29) 



for X » 1 jAbramowitz & Stegunlll972h . The resulting approxima- 
tion for the Wigner-3 j symbol overestimates the true value on av- 
erage by ^ 0.4% in the relevant ^-range. Weak lensing convergence 
bispectra resulting from this procedure are shown in Fig.[5]for the 
equilateral configuration. As in the case of the convergence spec- 
tra C,t(f), the bispectra Bi^{Cp) attain higher values in models with 
decaying CDM compared to those with stable CDM. This increase 
amounts to ^ 10%, and is observed independent of the equation of 
state of dark energy, where one observes a change of a few percent 
in the models with a cosmological constant relative to those with a 
varying equation of state. 

Lensing convergence bispectra for isosceles triangles are 
given in Fig.|6l for { = 10^. The plot suggests that the increase in 
power observed in decaying CDM models is present for all config- 
urations, if the background galaxy distribution has a high average 
redshift, which is in fact expected for a change of cosmology on 
the homogeneous level. Again, typical differences between mod- 
els with stable CDM and those with F = 1/3 amount to ^ 10%, 
irrespective of the opening angle of the triangle. 

Fig. |7] illustrates the configuration dependence of the weak 
lensing bispectra, for = 1000. In particular, the plot shows the 
configuration dependence variable Rc^{£\,{2), 



BiliJiJi) 



(30) 



for ^3 = 10^. Differences in cosmology are cancelled in first order 
in this expression, and for that reason the plot only shows the bis- 
pectrum configuration dependence for the fiducial ACDM model. 
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Figure 6. Isosceles convergence bispectra Bi^{C, cos ijj): ACDM (solid line), 
ArCDM with F = i (dashed line), ^CDM (dash-dotted line), and (^rCDM 
with r = ^ (dotted line), for the entire galaxy sample as a function of the 
cosine of the opening angle ij/. The angular scale is fixed to t{ = C2 = 10^. 
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Table 2. Characteristics of the weak lensing survey of the Dark UNiverse 
Explorer (DUNE), in terms of galaxy number n per sterad (corresponding 
to a galaxy density of 40/arcmin'), solid angle Af2, sky fraction /^j^y, survey 
depth 7,0 and intrinsic ellipticity dispersion cr^ of the galaxy sample. 



The tomography bins are chosen to contain equal numbers of back- 
;round galaxies, w hich are distributed in redshift z according to 
Small etal .11 1995h . 



p(z)Az = po 



— exp 

Zo 



-(-I |dz 

:o. 



wtth l = ^rf^l 



(34) 



where the distribution is characterised by the parameters /3 = | and 
Zo = 0.64 (corresponding to a median redshift of Zmsxi = 0.9, c.f. Ap- 
pendix |B](. We will consider 2-bin tomography, because any finer 
subdivision does not necessarily improve the bounds on cosmolog- 
ical parameters due to the high covariance of the signal originat- 
ing from diffe rent tomography bins, as discussed in iTakada & JainI 
( l2003bl.l2004l) . 




multipole order 



multipole order £2 



Figure 7. Convergence bispectrum configuration dependence {ii , ^2) for 
^3 = 10', without subdivision into tomography bins. Non-admissible trian- 
gle configurations violating the inequality |f, - < 4- ^ ti + l^j would 
occupy the left empty side of the plot. 



4 PARAMETER CONSTRAINTS 

We use a Fisher-matrix approach jTegmark et alj 1997h for estimat- 
ing the accuracy of the determination the cosmological parameters 
Q„,, cTg, the dark energy properties wq, and the CDM decay rate 
r for the weak lensing survey to be carried out by the DUNE experi- 
ment. The most important characteristics of DUNE are summarised 
in Table|2] 

4.1 Bispectrum covariances 

The observed bispectra -S(,yQ(/'i, ^2, ^3) are unbiased estimates of 
the true bispectra S(,jt)(^i, ^2> ^3)> 



B(ijk){il,(2,(i) - 5(;7«(^1> ^2> ^3)^ 



(35) 



because the intrinsic ellipticity distribution p(e)de of the back- 
ground galaxies is assumed to be skewless. But p(e)de has a fi- 
nite width cr^ , which im pacts on the observed power spectra Cij(£) 
( iKaiseill 9981 : Ihi]! 19991) . 



C(ij)(i) = Ciij){l) + 6ij — . 



(36) 



3.3 Bispectrum tomography 

In the case of bispectrum tomography, where the background galax- 
ies are divided in to two or more redshift bins (i,j,k), eqn. ([28} 
generalises to 



r 

fi(,7«(^i,^2,^3)= dx 
Jo 



D'i(x)B,(ki,k2,ki),0l) 



with indidual weighting functions 



(32) 



,X>Xi+i^ 

where the lensing-kemel weighted redshift distribution is given by 



dx p{z)— — - 

JmMx.X,) X 



(33) 



rij is the number of galaxies per steradian in the tomography bin 

(37) 



iii = n dx P(z)— , 



with the total number of galaxies per steradian n. The bispectrum 
covariance for to mography, which is diagonal in £„, p = 1, 2, 3, is 
approximated by ( [tMl200(]|:lTakada & Jairj|2003bll2004l) 



COV fl(,/H(C)6(;„„)(^„) 



A(A,^2,^3)~ 

■C'(;7)(^l)C(jm)(^2)C(*:n)(^3)> (38) 



7 .sky 



where the function A{ip) counts the multiplicity of triangle config- 
urations and is defined as 



MeiJ2J3) 



6, f 1 = ^2 — ^3 

2, ii = (j for / ^ j 
1, 1 12 + 



(39) 
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/sky denotes the fraction of the observed sky. It is worth noting that 
the bispectra 6(i ,^2,^3) are invariant under permutations of the 
bin indices. Because of that, there are 4 independent bispectra for 2- 
bin tomography, which can be conveniently indexed with the num- 
ber^ = 2' + 2' + 2*. 



4.2 Fisher matrices 

The Fisher matrix for the parameter space consisting of x 6 
{Sim, Ts, Wo, Wa, T, h] is constructed with ACDM as the fiducial cos- 
mological model. Particularly, for measurements of the weak lens- 
ing bispectrum S(,jt)(f|,^2,^3), one obtains: 



(^L,OVLB(,jj.)Z*(;„„)l j 



(Imii) 



a.m,,,) 



dxy 



(40) 



with (p e [li,{2,ti]. The summation is carried out with the condi- 
tion Ci < ^2 ^ ^3, such that every triangle configuration is counted 
once, between the limits {^^^ = 10^ and ^,„ax = 2 x 10^. The indices 
J, k and /, m, n run over the tomography bin s. 

Following Imada & JainI j2003bl . l2004h . we use binned sum- 
mations in the multipoles li and £2, but carry out an unbinned sum- 
mation in ti in order to account for the vanishing Wigner-3 j symbol 
if Yip ip is an odd number and for the sign change of the Wigner-3 j 
symbol depending on whether Y,p ^p mod 4 vanishes or not. Using 
A/'i = A/'2 = 25, the computational load is reduced by three or- 
ders of magnitude, as we have to compute 2x10' triangles instead 
of 3 X 10'" triangles, for l^^^ = 2 x 10''. Using binned summation 
in two variables is well justified, given the smooth variation of the 
bispectrum illustrated in Fig. [T] The scale on which the numerical 
derivatives are computed corresponds to a 5% variation for erg 
and h and to a 10% variation for the dark energy parameters wo, Wa 
and r. 

We would like to emphasise at this point that our Fisher- 
analysis implicitly assumes priors on (i) spatial flatness, n„, + Q.^ = 
1, which affects all geometrical measures, (ii) the baryon density 
Q.t, which adds a correction to the shape parameter, and finally (iii) 
the slope n, of the CDM spectrum P{k) for small wave numbers 
k. We neglect the influence of Q.t, because the baryons cause only 
a minor correction to the shape parameter, and keep = 1 and 
Q.,„ + = 1 fixed, as they are both generic predictions of inflation 
and well tested by CMB observations. The parameter space consid- 
ered here is motivated by the fact that n,„, w'o, and F all increase 
the lensing signal, either by their influence on the growth equation 
or by appearing in the Poisson equation, and are naturally degen- 
erate with (Tg. As additional priors to the lensing measurement, we 
assume Icr-errors from PLANCK CMB observatio ns with the nu - 
merical values Ish = 0.13 and Acrg = 0.01 (Eisenst ein et alJI 19991 ). 
The priors are assumed to be diagonal, f ™^ = Sf^y/cr'^, and added 
to the lensing Fisher matrix F^^, 



(41) 



because the likelihoods of independent measurements can be mul- 
tiplied and their ;j'^-functions added. 

The diagonal elements of the inverse Fisher matrix give the 
Cramer-Rao bound on individual parameters. 



(42) 



which are compiled in Table [3] for 2-bin weak lensing tomogra- 
phy, with the enhancement of the measurement by including CMB- 
priors on erg and on h. Quite generally, weak lensing bispectrum 
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Table 3. Expected accuracy on the cosmological parameters from 2-bin 
weak lensing tomography, with ACDM as the fiducial cosmology. The sec- 
ond column adds a CMB prior on o"g and h to the Fisher-matrix. 



tomography with DUNE in conjunction with CMB priors provides 
percent errors on Q.,,,, erg and h, errors of the order of 20% on w'o 
and F, but a weak constraint o n Wg, which is due to the fact that 
in comparing to the analysis bv lTakada & Jairj ( l2004h . the DUNE 
galaxy sample has a much lower median redshift (0.9 compared to 
1.5) and hence a weaker lever arm on w^. 

The constraint on AF ^ 0.13 translates to a lower bound on 
the CDM lifetime of more than 7.7 Hubble times, corresponding 
to tj- i 75.3//1 Gyr. Limits on the CDM particle lifetime from lens- 
ing measurements of that order of magnitude may serve to exclude 
certain particle candidates, but we should emphasise that we inves- 
tigate only a single decay channel with an unknown branching ratio 
b: If CDM has two decay modes with the probabilities b for decay- 
ing into dark energy and 1 - fo for decaying into other products, with 
corresponding decay rates F^ (constrained by lensing) and Fx, the 
total decay rate is given by F = bT^ + (1 - b)Tx- 

The^--value for pairs of parameters (x^, Xy) can be computed 
from the inverse (F"')^,, of the Fisher matrix jMatsubara & Szalavl 
|2002|) . 



Axp 

AXy 



{F-%p (F-%, 



Axp 

AXy 



where Ax^ 



(F-X (F-\ 
- x^'^^'^. Specifically, the contour at 



(43) 



A;r' = -21n erfc|-^j, (44) 

encloses a fraction of erf(n/ V2) of parameter space, which corre- 
sponds to a confidence level of ncr. erf(x) and erfc(x) = I - erf(x) 
are the err or function and the complem entary error function, re- 
spectively jAbramowitz & StegunI [l972[) . The correlation coeffi- 
cient r,,y is defined as 



r.y 



{F-')py 



yjiF-%,{F-')r 



(45) 



and describes the degree of dependence between the parameters x,, 
and Xy by assuming numerical values close to for independent, 
and close to I for strongly dependent parameters. Constraints on 
pairs (Xp,Xy) of cosmological parameters from the Fisher- analysis 
are compiled in Fig. [8] along with the respective correlation coeffi- 
cient r-fiy. 

The CDM decay constant F is negatively correlated with Q.„, 
because a simultaneous decrease of the matter density today as well 
as a increase in the decay constant leads to the same averaged mat- 
ter density for the lensing signal. The degeneracies with wq and w„ 
are such that wo has to be de- and Wa increased in order to have the 
same lensing signal in models with non-vanishing F. The positive 
correlation of F with trg is a projection effect in the marginalisation 
driven by the tight constraints between Q„,, wq and h, and the classi- 
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Figure 9. Fisher matrix contraints on the triplet (erg, wq, F), from weak lens- 
ing tomography including a CMB-prior on erg and h. The ellipsoids coiTe- 
spond to lo" and 2tT intervals. 



Figure 10. Fisher matrix contraints on the triplet (n„,,(T8,F), from weak 
lensing tomography including a CMB-prior on erg and h. The ellipsoids 
correspond to Itr and 2o" intervals. 



cal degeneracies between Q.,,, and erg (the power spectrum normali- 
sation and the strength of the gravitational potentials) and between 
n„, and h (the definition of the shape parameter) are recovered. 

We would like to point out that the effective equation of state 
parameter w^ff may cross the boundary towards phantom models, 
for which w < -1, 



Weff = W 



3h p0 



Weff 



3 



(46) 



with H{a) = HQh(a), if decay is considered. Therefore, we mark 
the forbidden regions of parameter space in the relevant plots. 
Generalising eqn. ( 143 l l to three parameters (x^, Xy, Xp) yields 
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(47) 



In a 3-dimentional parameter space, the contours at Axj - 2.6124 
and Axl - 4.2222 correspond to significance levels of Icr and 2cr, 
respectively. Figs.l9landllOlsummarise simultaneous constraints on 
the triplets (erg, wq, T) and (Q.,,,, (Tg, F), respectively. The increase in 
lensing signal by either increasing trg or the dark energy equation 
of state H'o is illustrated in Fig. |9l with the weak degeneracy that 
models with decaying CDM require higher values of erg, for lens- 
ing at low redshift. Fig.fTO] on the contrary, shows that erg can be 
decreased in models with decay if fi„, is increased at the same time. 



5 SUMMARY 

In this paper, we investigate the capability of the DUNE experiment 
to constrain the decay of dark matter into dark energy from the 
observation of the weak convergence bispectrum, in a tomographic 
measurement. 

(i) Our cosmological model with CDM decaying into dark en- 
ergy may provide an alternative explanation of the coincidence 
problem, i.e. the fact that the cosmic expansion is dominated by 
dark energy after structure formation. This is achieved by choos- 
ing a small value of the decay constant F. Decaying CDM naturally 
influences the Hubble function by providing dark energy, and by in- 
troducing a faster scaling of the matter density compared to p oc a^^ 
in models with stable CDM. 



(ii) We have neglected the contribution of baryons, whose den- 
sity would just decrease oc tr\ because they are stable particles. 
This should not be a serious limitation, however, due to the fact 
that the baryon fraction = Qi/O,,, =^ 0.16 has a small numerical 
value. Models with decaying CDM have the interesting property, 
however, that the baryon fraction = ptlPm is monotonically in- 
creasing with cosmic time. 

(iii) The growth of structure is influenced in two ways: Firstly by 
the non-standard scaling of the density parameter Q.„, which causes 
stronger gravitational potentials in the past in models with decay- 
ing CDM, and secondly by affecting the logarithmic derivative of 
the Hubble function, which is smaller in decaying models. These 
changes are naturally degenerate with models for dark energy with 
a varying equation of state. We work in the limit of the dark en- 
ergy sound speed being close to c, such that dark energy can be 
considered homogeneous even though it originates by decay from 
clustered dark matter. 

(iv) Similarly to the growth of structure, the coupling to of light 
to the matter in gravitational lensing is affected by the higher value 
of f2,„, which causes the gravitational potential O induced by the 
overdensity field S to be stronger in decaying models. With this 
mechanism, models with decaying CDM can provide an explana- 
tion for the high values of trg required by cosmic shear experiments, 
and can be reconciled with measurements of erg from CMB obser- 
vations. 

(v) Weak lensing bispectrum tomography yields a relative ac- 
curacy on the determination of the dark energy equation of state 
parameters wo, vvo, the CDM decay rate F of a tenth and percent 
accuracy on the parameters 0,„, erg, if CMB priors are included. 
The CDM decay rate F is naturally degenerate with the equation of 
state parameters wo and Wa- The Fisher-analysis provides an upper 
bound on F, or equivalently, a lower bound on the CDM lifetime, 
fr = 1/F > 1.1 1 Ha ^ 75.3 Gyi/h. This limit might be useful to 
exclude certain CDM particle candidates, although we should em- 
phasise here that we investigate a single decay channel only, and 
that the branching ratio of this particular channel would need to be 
known as well. 

(vi) Comparing the lensing constraints on F to supernova con- 
straints shows that differences in luminosity distance small: There 
is a 0.01% difference between ApCDM and ACDM and a 1% dif- 
ference between OpCDM and ACDM at z = 0.1. Moving to higher 
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Figure 8. Constraints on cosmological parameters and the CDM decay rate from the Fisher-matrix analysis from weak lensing tomography. The choice for 
the fiducial cosmological model is £1,,, = 0.25, erg = 0.8, wq = -I, Wa = 0,T = and h = 0.72 (marked by a circle). The ellipses correspond to lo" . . . 4cr. The 
correlation coefficient r^,, of the pai'ameter pair (jc^ , JCv) is given in the upper right comer of each panel (to two digits), and the the excluded region of parameter 
space due to the phantom constraint eqn. )46t is below the line in the wq-V plot. 



redshifts confirms the trend that the equation of state has a stronger 
influence on the luminosity distance than the CDM decay rate. At 
z = 1, one observes a 1% difference between ApCDM and ACDM, 
but 7% difference between the stable dark energy models <I)rCDM 
and ACDM. 

(vii) It is worth noting that models with decaying CDM are 
genuinely different from dark energy models concerning structure 
growth and observations which use gravitational interaction such 
as gravitational lensing. While it is always possible to construct an 
equation of state w{a) for a dark energy model with stable CDM 
that gives the identical Hubble function as a model with decaying 
CDM, the evolution of the density parameter Cl,„(a) and the growth 
function D+(a) breaks this degeneracy which would be different in 
the two cases. For that reason, the combination of probes of cosmic 
structure growth and the expansion history are able to distinguish 
between the two families of models, which would not be possible 
with e.g. supernova observations alone. 

In addition, we plan to provide prospective constraints from 
the integrated Sachs-Wolfe effect, which is promising as it orig- 
inates at higher redshifts, and directly measures the derivative 
dlnp„,(a)/dlna, which is t -3 in decaying CDM models. In sum- 
mary we would like to stress that CDM-decay would be an elegant 
solution to the coincidence problem, and that its central parameter 
r is well measurable by future lensing surveys. 
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Figure Al. Comparison of the dissipation term Q{a) = 3 + dlnff/dlna 
(thick lines) and the source term S(a) = Q„,{a) (thin lines) in the growth 
equation, for ACDM (solid line), ApCDM with F = ^ (dashed line), 0CDM 



(dash-dotted Une), and 0rCDM with T = j (dotted line). 
Wang L., Steinhardt P J., 1998, ApJ, 508, 483 

APPENDIX A: GROWTH FUNCTION 

The growth function D+(a) is a solution to the differential equation 

-n„(a)0+(a). (AI) 



d^ 1 / dln/f 

— D+(fl) + - 3 + — 

aa-^ a \ din a 



d 3 



da ' 2a^ 

Fig. |AT] compares the source term 5(a) = Q.„Xa) with the dissipa- 
tion term Q(a) = 2- g = 3-l-dln H/d In a for the exemplary cos- 
mologies. Dark energy with a variable equation of state increases 
the dissipation term and decreases the source term of the growth 
equation, and suppresses growth relative to ACDM. CDM decay 
is able to compensate this effect, by enhancing growth due to de- 
creasing the dissipation term Q(a) even below the value of 1.5, and 
a slower increase at late times. At the same time, the matter density 
remains closer to the critical density for a longer period of time, 
until the CDM decay causes it to decrease to the current value. 

The near-degeneracy of the growth functions D+(a) in ACDM 
and 0rCDM are explained by the fact that they show very similary 
evolutions of Cl„,(a), and up to z = 1, the damping term Q(a) is 
almost identical between the two models. At smaller redshifts, the 
damping as well as the source term are both smaller compared to 
ACDM, with a compensating effect on the evolution of D+(a). 



APPENDIX B: REDSHIFT DISTRIBUTION 

The unit-normalised redshift distribution p{z)dz dSmail et al.ll995b . 



p(z)dz = po\ — \ exp 



£[ dz wtth l = ^r 3 



the cumulative distribution P(z), 
P(z)- 



f dz'pizl 
Jo 



(Bl) 



(B2) 



and the redshift boundaries for 2-bin (z = 0.90) and 3-bin (z = 0.72 
and I. II) tomography are shown in Fig. IB II The parameters are 
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Figure Bl. Differential distribution p(7)Az (thick solid line) and cumulative 
distribution P{z) = dz'piz') (thick dashed line) of the lensed galaxies 
observed by DUNE, with redshift bins for 2-bin (thin solid line) and 3-bin 
(thin dashed line) weak lensing tomography. 

Zo = 0.64 (corresponding to a median redshift of Zmed = 0.9) and 
j8 = I . These bins are chosen in such a way that there are equal 
number of background galaxies in each tomography bin. This has 
the consequence that the noise contribution cr^/n; to the spectral 
measurements from each bin is equal. Although convenient, this 
choice of redshift binning is by no means optimised for providing 
the maximal accuracy in the measurement of cosmological param- 
eters. It should be possible, however, to perform an iterative proce- 
dure with a coarse determination to begin with, and a subsequent 
choice of tomography bins which maximises the accuracy. 
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